First-principles Calculation of Superconductivity in Hole-doped LiBC: Tc = 65 K 
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The lattice dynamical properties of Li^BC are calculated for several values of x using density 
functional perturbation theory. We find that the electron-phonon coupling parameter A increases 
P5 . monotonically with decreasing a; to a maximum value of 1.4 for x — 0.125 owing to the increasing 

O ' radius of multiply-nested Fermi surface cylinders. The B-C bond-stretching phonon modes have 

O ' frequencies which are 28% higher than the equivalent modes in MgB2. This combination results in 

5^ , a Tc of about 65 K for x = 0.5. 

Oh' 

c/3 ■ It has been suggested recently that hole doped LiBC may have a superconducting transition temperature similar 
(_i [ to or exceeding that of the isovalent compound MgB2 |jl|, g]. Experimental confirmation of this would be of great 
^2 ' significance since MgB2 already has an exceptionally high transition temperature of ^ 40 K. Since the discovery of 
superconductivity in MgB2 |^ there has been a flurry of interest in it and compounds with similar features, namely, 
layered crystals with weak interlayer interactions, a corresponding high density of states at the Fermi energy and 
strong in-plane covalent bonds which are partially filled |, |, ||, |, |, |, |Il|, |l2|, ^ ^ [l6[ LiBC has 
a hexagonal structure with alternating lithium and boron-carbon planes. Although it is semi-conducting owing to 
Q ' the alternating B-C atoms within each plane [Q, the addition of holes, likely to be accomplished by the removal 
' of Li atoms, cause it to become metallic and thus potentially superconducting. To date no such effect has been 
^ • observed experimentally which is the motivation for this article in which we report a first-principles calculation of the 
^ ' electron-phonon coupling in Lij;BC. 

To obtain accurate results we employ a plane wave method with the analytic pseudopotentials of Hartwigsen, 
Goedecker and Hutter using the parameters quoted in their paper jl^ . A kinetic energy cut-off criterion of 24 Hartree 
i was placed on the plane waves to ensure good convergence in terms of basis size. The density functional linear 
' response method of Giannozzi et al. p9| ] was used to calculate thephonon spectra and corresponding self-consistent 
linearized potentials within the local density approximation |2^, |2^. Dynamical matrices and linearized potentials 
were calculated in the irreducible Brillouin zone on a 4 x 4 x 2 g-point mesh using the ground-state wave functions 
obtained on a 8 x 8 x 4 fc-point mesh. To incorporate hole doping we removed the appropriate number of electrons 
"j^ ■ from the unit cell and added a positive uniform background to maintain charge neutrality. This approximation is 
I justified from the results of Rosner, Kitaigorodsky and Pickett (RKP) |l[ who noted an almost rigid shift in the Fermi 
I . level with little rcnormalization of the electronic density of states upon removal of lithium atoms. The phonon spectra 
■ and electron-phonon coupling results were calculated for x — 0.125, 0.25, 0.375, 0.5, 0.625, 0.75 and 0.875. We also 
Ch , performed a ground-state calculation on an even denser 16 x 16 x 8 fc-point mesh for producing the electron-phonon 
^ • coupling matrices as well the Fermi surface plots in Fig. 0. 

\ In Fig. ^ we plot the calculated phonon dispersions along high symmetry lines for each doping level. As is evidenced 
^ from the lack of singularities in the phonon curves, the material is found to be stable for all doping concentrations 
k> [ within the approximation of a uniform background. The most noteworthy feature in the spectra is the abrupt dip 
. in the topmost branches between the F — M, A' — F and A ~ L points. As with MgB2, these branches have the E2g 
d ' symmetry character and the dip arises from the exceptionally strong nesting between concentric cylinders | ]l7[ of the 
Fermi surfaces which are oriented in the direction of the c-axis as shown in Fig. |l]. These nesting modes correspond 
to in-plane bond stretching between the boron and carbon atoms. With decreasing x (i.e. increasing doping) the 
sharpness of the dip becomes less pronounced although the midpoint moves away from the central axis. Also of 
interest is the level attraction between certain phonon branches for x < 0.5. This is reflected in the Fermi surface 
cylinders which change from being nearly circular for low doping to hexagonal for high doping. This flattening of 
the sheets results in additional translational symmetries in the effective potential. The C-B bond-stretching modes 
have phonon frequencies around 1000 cm^^ which is considerably higher than the B-B frequencies of MgB2 which are 
about 725 cm^^ 111. 
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The phonon density of states are plotted in Fig. ^. At low doping levels there is a double peak structure between 
200 and 400 cm^^ as well as peaks at 620 and 1000 cm~^. As the doping level is raised the lower peaks merge together 
and the higher frequency peaks become much less pronounced. The gap-like structure in the density of states between 
400 and 600 cm^^ disappears with increased doping. 

We now turn to calculation of the superconducting transition temperature. Details of Migdal-Eliashbcrg theory 
which is the basis for these calculations can be found in Rcf. The electron-phonon interaction matrix is 

determined from the linearized self-consistent potential (5Vq by 



0k^y, (1) 
where Mi is the mass of the i'th ion, is a single-particle state and e^j-.i is the eigenvector of phonon mode qj and 
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FIG. 3: Calculated phonon density of states for various hole dopings. 
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FIG. 4: Calculated electron-phonon spectral function for various hole dopings. 



atom i. From this we obtain the phonon linewidth arising from the electron-phonon interaction by averaging over the 
Fermi surface 

7qj- = 27rujqj \9k+qv',kv\^^i^i^" " eF)(5(ek+qi,' - ep), (2) 

where is the Kohn-Sham eigenvalue and e-p is the Fermi energy. The central quantity in superconductivity 
theory is the electron-phonon spectral function a'^F{uj). This can be written in terms of the phonon linewidth by 

a'Fic) = -i- J2 —Si^ - (3) 

where Np is the density of states at the Fermi energy. The total electron-phonon coupling parameter is defined as 
the first reciprocal moment of a^F{u!) 

X = 2r^d.. (4) 







Finally Tc can be estimated from the McMillan equation |24 



1.2 ^ V-^-M*(l + 0.62A) 

with fi* , the Morel- Anderson pseudopotential, representing Coulombic repulsion and the logarithmic average frequency 
given by 

(2 f°° ^ , Anuj] 
u;iog = exp|-y aF{uj) — j. (6) 

The three main components in determining Tc are therefore the density of states at the Fermi energy Np, the 
logarithmic average frequency wiog and the strength of the electron-phonon coupling on the Fermi surface. By 
plotting a^F^Lu) we can gauge the relative strength of the coupling. This is done in Fig. ^. For most phonon- 
mediated superconductors the phonon density of states and a'^F{u}) differ only slightly j2^. This is not the case here: 
the modes above 600 cm"-'^ which dominate in a'^F^uj) are suppressed in the phonon density of states in Fig. ^ This 
is particularly evident for the peak at x — 0.5 and uj « 600 cm~^ which corresponds to the weakly dispersed i?2g 
branch in the F — ^ direction. We plot in Fig. || N-p, wiog, A, and Tc using an empirical value for /i* of 0.09 which 
was found to reproduce Tc in MgB2 The density of states at the Fermi energy increases sharply from x — 1 to 
X = 0.75. Beyond this N-p increases more gradually with a slight upturn at x = 0.5 and has values between 15 and 
21 states/spin/Hartree/unit cell. As hole concentration is raised w\og drops by about 25% from 802 K at x = 0.875 
to 608 K at a; = 0.125. There is also an inflection point a,t x = 0.5 in this curve. The total electron-phonon coupling 
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FIG. 5: Properties calculated as a function of Li concentration x: density of states at Fermi energy in units of 
states/spin/Hartree/unit cell (a), logarithmic average frequency (b), electron-phonon coupling parameter A (c), and supercon- 
ducting transition temperature from the McMillan formula (d). The circles represent calculated points and the interpolating 
lines are cubic splines. 

parameter A increases over the whole range to a maximum value of 1.43. This is a direct result of the increase with 
1 — a: of the Fermi surface nesting vector 2fcF which maps the opposing walls of the cylindrical sheets shown in Fig. 
|. We again note an inflection point in A for x = 0.5. From proportionality arguments, RKP obtained A « 1.5 for 
X = 0.5 which is larger than the 1.20 found here. Nonetheless, this is substantially larger than A = 0.82 found in 
MgBa H H, Despite the fairly large increase in A over the whole range of x, Tc reaches a plateau of 65 K at 
X = 0.5. The reason for this is the drop in wiog as a function of 1 — a; which enters as a proportionality in Eq. ^ as 
well as the saturation of iVp for x < 0.75. The highest value for Tc is 68 K at a; = 0.125, this is substantially higher 
than that observed for MgB2 but lower than the ~ 100 K estimated by RKP. 

In conclusion, we find from our first-principles calculations that hole-doped LiBC is a superconductor with ex- 
ceptionally strong in-plane coupling. As with MgB2 this is due to a small number of phonon modes close to the 
r point which scatter electrons at the Fermi surface from one side of a cylindrical sheet to another. Although the 
superconducting transition temperature is much higher than in MgB2, the decrease in the average phonon frequency 
as the doping concentration is raised limits the potential gain in Tc from the increase in A and we note that hole 
doping beyond x = 0.5 will not result in any substantial improvement Tc- Nonetheless, going by these results, Li^^BC 
may be an impressive superconductor well worth further investigation. 
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